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Abstract 

We study the branching ratios, the direct CP asymmetries in B ^ /^{*)^{*) decays 
and the polarization fractions of — > K*K* decays by employing the QCD factorization in 
the minimal supersymmetric standard model with R-parity violation. We derive the new 
upper bounds on the relevant R-parity violating couplings from the latest experimental 
data of i? ^ K^*^ K^*\ and some of these constraints are stronger than the existing 
bounds. Using the constrained parameter spaces, we predict the R-parity violating effects 
on the other quantities \n B ^ decays which have not been measured yet. 

We find that the R-parity violating effects on the branching ratios and the direct CP 
asymmetries could be large, nevertheless their effects on the longitudinal polarizations of 
B — > K*K* decays are small. Near future experiments can test these predictions and 
shrink the parameter spaces. 
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1 Introduction 



The study of exclusive hadronic S-meson decays can provide not only an interesting avenue to 
understand the CP violation and flavor mixing of the quark sector in the standard model (SM), 
but also powerful means to probe different new physics (NP) scenarios beyond the SM. Recent 
experimental measurements have shown that some B decays to two light mesons deviated from 
the SM expectations, for example, the tttt, ttK puzzle P and the polarization puzzle in B VV 
decays [2]. Although these measurements represent quite a challenge for theory, the SM is in no 
way ruled out yet since there are many theoretical uncertainties in low energy QCD. However, 
it will be under considerable strain if the experimental data persist for a long time. 

Among those NP models that survived electroweak (EW) data, one of the most respectable 
options is the R-parity violating (RPV) supersymmetry (SUSY). The possible appearance of 
the RPV couplings [Hj, which will violate the lepton and baryon number conservation, has 
gained full attention in searching for SUSY [HE]. The effect of the RPV SUSY on B decays 
have been extensively investigated previously in the literatures IZj ? and it has been proposed 
as a possible resolution to the polarization puzzle and the vrvr, nK puzzle |H]. The pure penguin 
B K^*^K^*^ decays are closely related with the puzzles which are inconsistent with the SM 
predictions, and therefore are very important for understanding the dynamics of nonleptonic 
two-body B decays, which have been studied in Refs. jHj. If the RPV SUSY is the right model 
to resolve these puzzles, the same type of NP will affect B k(*)k^*^ decays. In this work, we 
shall study the RPV SUSY effects in the B ^ K^*)Ki*) decays by using the QCD factorization 
(QCDF) approach jTUj for hadronic dynamics. The B }{(*)}{(*) decays are all induced at 
the quark level by 6 ^ dss process, and they involve the same set of RPV coupling constants. 
Using the latest experimental data and the theoretical parameters, we obtain the new upper 
limits on the relevant RPV couplings. Then we use the constrained regions of parameters to 
examine the RPV effects on observations in the B — > K^*^K^*^ decays which have not been 
measured yet. 

The paper is arranged as follows. In Sec. 2, we calculate the CP averaged branching ratios, 
the direct CP asymmetries oi B ^ and the polarization fractions in i? K*K* 

decays, taking account of the RPV effects with the QCDF approach. In Sec. 3, we tabulate 
the theoretical inputs in our numerical analysis. Section 4 deals with the numerical results. 
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We display the constrained parameter spaces which satisfy all the experimental data, and then 
we use the constrained parameter spaces to predict the RPV effects on the other observable 
quantities, which have not been measured yet in B system. Section 5 contains our 

summary and conclusion. 



2 The theoretical frame for B K^*^ K^*^ decays 
2.1 The decay amplitudes in the SM 

In the SM, the low energy effective Hamiltonian for the AB = 1 transition at the scale n is 
given by [TT] 



r 10 ^ 

'^e// ~ ~~/^ X/ ^P] ^^Ql + ^2^2 + \CiQi + Cj^Qj^ + CsqQsq >+h.C. 
V ^ p=u,c [ 1=3 ) 



(1) 



here \p = VphV*^ for 6 — > g transition (p G {m, c}, q G {rf, s}) and the detailed definition of the 
operator base can be found in [TT] . 

Using the weak effective Hamiltonian given by Eq.Q, we can now write the decay ampli- 
tudes for the general two-body hadronic B — > M1M2 decays as 

A^^{B ^ M1M2) = {MiM2\H^fj\B) 

= %EE^C'.(/i)(MiM2|Q,(/x)|E). (2) 

The essential theoretical difficulty for obtaining the decay amplitude arises from the evaluation 
of hadronic matrix elements {MiM2\Qi{^)\B) . There are at least three approaches with different 
considerations to tackle the said difficulty: the naive factorization (NF) ^21^], the perturbative 
QCD [m, and the QCDF The QCDF developed by Beneke, Buchalla, Neubert and 

Sachrajda is a powerful framework for studying charmless B decays. We will employ the 
QCDF approach in this paper. 

The QCDF ^01 allows us to compute the nonfactorizable corrections to the hadronic matrix 
elements {MiM2\0i\B) in the heavy quark limit. The decay amplitude has the form 



(3) 
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Figure 1: The next to leading order nonfactorizable contributions to the coefficients a^. 




Figure 2: The weak annihilation contributions to the coefficients 6f . 



where the effective parameters af including nonfactorizable corrections at order of a^. They are 
calculated from the vertex corrections, the hard spectator scattering, and the QCD penguin 
contributions, which are shown in FigH The parameters 6f are calculated from the weak 
annihilation contributions as shown in Fig|21 

Under the naive factorization (NF) approach, the factorized matrix element is given by 



A 



M1M2 



^ (M2|fe7/.(1 - l,)q-m{M,\{h^{l - 75)gi)|5). 



(4) 



In term of the decay constant and form factors [T3], AmiM2 ^-^e expressed as 
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2 T7B-*Mi, 
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M2/' 



-ifM^mM^ 



2 AB~*Ah, 
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if Ml = P, M2 = P, 
if Ml = P, M2 = V, 
if Ml = 1/, M2 = P, 



{el ■ sDimn + m^JAf ^""Hmi.J - (el ■ pB)ie*, ■ pb) 



2A. 



+^e^lua/3£2^l PbPi 



2^ 



if Ml = V, M2 = V, 



(5) 



where P(V) denote a pseudoscalar (vector) meson, Pb{^b) is the four-momentum(mass) of the 
B meson, 771^. is the masses of the Mj mesons, and e* is the polarization vector of the vector 
mesons M,-. 



Following Beneke and Neubert (TH], coefficients af can be split into two parts: af = of j+af 
The ffist part contains the NF contribution and the sum of nonfactorizable vertex and penguin 
corrections, while the second one arises from the hard spectator scattering. The coefficients 
read [15] 
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where = as(//), Cp = {N^ — l)/{2Nc), Nc = 3 is the number of colors, and Nm2 = 1(0) for 
M2 is a pseudoscalar (vector) meson. The quantities Vm2, -f^A/iA/2) P^h 25 ^Ah 3' PhifY PmT 
consist of convolutions of hard-scattering kernels with meson distribution amplitudes. Specifi- 
cally, the terms Vm2 come from the vertex corrections in Figma)-md), -Pm2,2 -^A/2,3 (P^h^ 
and Pm^Y) arise from QCD (EW) penguin contractions and the contributions from the dipole 
operators as depicted by FigHl^e) andlH^f). HmiM2 is due to the hard spectator scattering as 
Fig^g) and^h). For the penguin terms, the subscript 2 and 3 indicate the twist 2 and 3 
distribution amplitudes of light mesons, respectively. Explicit forms for these quantities are 

relegated to Appendix A. 

We use the convention that Mi contains an antiquark from the weak vertex, for non-singlet 
annihilation M2 then contains a quark from the weak vertex. The parameters 6f = 6f (Mi, M2) 
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in Eq.Q correspond to the weak annihilation contributions and are given as ^7] 



&i(Afi, M2) = -^Ci^KMi, M2), ?>2(Mi, A/2) = -^C2Ai(Mi, M2), 

&f (Ml, M2) = [Cs^KMi, M2) + C5 (^^(Mi, M2) + 4(Mi, M2)) + NcC^Ai{Ah,M2) 
6?(Afi, M2) = [CMNh,M2) + C6^^(Mi, A/2)] , 

6f'^'"(Afi, Af2) = [CgAUAfi, A/2) + C7 (a^(A/i, Afs) + A^(Afi, Afs)) + iVcC84(A^i- ^2) 
65'^'"(Mi, A/2) = -ll^ [CioAl(A/i, A/2) + C8^^(A/i, M2)] , (7) 

the annihilation coefficients (61, 62), 64) and (63''^"', Tf^^^) correspond to the contributions of 
the tree, QCD penguins and EW penguins operators insertions, respectively. The explicit form 
for the building blocks A*/(Mi, M2) can be found in Appendix A. 

With the coefficients in Eq.® and ((Zj), we can obtain the decay amplitudes of the SM 

If 

part) and ^f^^ (the subscript "a" denotes the annihilation part). The SM part amplitudes of 
B K(*)Ki*) decays are given in Appendix B. 



part A^'^^ (the subscript "/" denotes the part without the contribution from the annihilation 



2.2 R-parity violating SUSY effects in the decays 

In the most general superpotential of the minimal supersymmetric Standard Model (MSSM), 
the RPV superpotential is given by 



where L and Q are the SU(2)-doublet lepton and quark superfields and E'^, Lf^ and Z)'^ are the 
singlet superfields, while i, j and k are generation indices and c denotes a charge conjugate field. 

The bilinear RPV superpotential terms fiiLiHu can be rotated away by suitable redefining 
the lepton and Higgs superfields ^Hl- However, the rotation will generate a soft SUSY breaking 
bilinear term which would affect our calculation through penguin level. However, the processes 
discussed in this paper could be induced by tree-level RPV couplings, so that we would neglect 
sub-leading RPV penguin contributions in this study. 

The A and A' couplings in Eq.® break the lepton number, while the A" couplings break 
the baryon number. There are 27 X[j,^ couplings, 9 Xijk and 9 X'-ji^ couplings. X[ij]k are antisym- 
metric with respect to their ffist two indices, and X'-y^ are antisymmetric with j and k. The 
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antisymmetry of the baryon number violating couplings A^^^j in the last two indices implies 
that there are no X'l^f. operator generating the b sss and h ddd transitions. 







yi \i* \ X A' \'* 

^ijk ^imn \ / ijk ^imn 

dk dn d„ 

Figure 3: Sleptons exchanging diagrams for nonleptonic B decays. 






A" A"* 



Figure 4: Squarks exchanging diagrams for nonleptonic B decays. 

From Eq.(jH)), we can obtain the following four fermion effective Hamiltonian due to the 
sleptons exchange as shown in FigOl 

id„^YPRdl)s{UkltiPLUj)8, 

{dr,^l^PRdlMdk-i^PLdj)8. (9) 
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The four fermion effective Hamiltonian due to the squarks exchanging as shown in FigH are 

kh - {ua^PRUj)8{dkj^PRdi)s 
{dkl^PRUj)i{uaf,PRdi)i - {dkYPRUj)8{ua^PRdi)s } , 
{da^PRdj)i{dk-i^,PRdi)i - ida^PRdj)sidk-f^PRdi)s\ ■ (10) 
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In Eqs.© and Pl 



1-7; 



1+75 



2 = -^^V = a^m'i,) ^'^'^ /5o = 11 - Inf. The subscript 



for the currents (j^)i,8 represents the current in the color singlet and octet, respectively. The 
coefficients r/~^/^o and ?^~^/^o are due to the running from the sfermion mass scale rrij (100 GeV 
assumed) down to the m?, scale. Since it is always assumed in phenomenology for numerical 
display that only one sfermion contributes at one time, we neglect the mixing between the 
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operators when we use the renormahzation group equation (RGE) to run Tt^jf down to the low 
scale. 

The RPV amplitude for the decays can be written as 

A^{B^MiM2) = (MiM2\7^}f\B) . (11) 

The product RPV couplings can in general be complex and their phases may induce new 
contribution to CP violation, which we write as 

here the RPV coupling constant A G {A, A', A"}, and is the RPV weak phase, which may 
be any value between — vr and vr. 

For simplicity we only consider the vertex corrections and the hard spectator scattering in 
the RPV decay amplitudes. We ignore the RPV penguin contributions, which are expected 
to be small even compared with the SM penguin amplitudes, this follows from the smallness 
of the relevant RPV couplings compared with the SM gauge couplings. The bounds on the 
RPV couplings are insensitive to the inclusion of the RPV penguins [20] • We also neglected the 
annihilation contributions in the RPV amplitudes. The R-parity violating part of the decay 
amplitudes A^^ can be found in Appendix C. 

2.3 The total decay amplitude 

With the QCDF, we can get the total decay amplitude 

A{B M1M2) = A^f^'{B M1M2) + Al^{B ^ M1M2) + A^^ {B ^ MM). (13) 

The expressions for the SM amplitude Aj^ and the RPV amplitude A^^ are presented in 
Appendices B and C, respectively. From the amplitude in Eq. |T!?jl . the branching ratio reads 

}3{B ^ M1M2) = \A{B ^ M1M2) f S, (14) 

where S" = 1/2 if Mi and M2 are identical, and S = 1 otherwise; tb is the B lifetime, |j9c| is the 
center of mass momentum in the center of mass frame of B meson, and given by 

\Pc\ = - i^Mi + "^Af2)^][?Ti| - {niMi - "iA/a)^]- (15) 
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The direct CP asymmetry is defined as 



B{B ^ /) - B{B ^ f) 



B{B^f) + B{B^fy ^^^^ 
In the B VV decay, the longitudinal polarization fraction is defined by 

-n 1/112,1/112,1/112' v-^ ' / 



r 

where Aq{A±) corresponding to the longitudinal(two transverse) polarization amplitude(s) for 
B ^VV decay. 

3 Input Parameters 

A. Wilson coefficients 

We use the next-to-leading Wilson coefficients calculated in the naive dimensional regularization 
(NDR) scheme at mf, scale 

Ci = 1.082, C2 = -0.185, C3 = 0.014, C4 = -0.035, 6*5 = 0.009, 

Cg = -0.041, Cj/ae = -0.002, Cg/tte = 0.054, Cg/a, = -1.292, 

Cio/tte = 0.263, 6*7;^^ = -0.299, = -0.143. (18) 

B. The CKM matrix element 

The magnitude of the CKM elements are taken from [21]: 

\Vud\ = 0.9738 ± 0.0005, |K.| = 0.2200 ± 0.0026, |Kb| = 0.00367 ± 0.00047, 

iKdl = -0.224 ±0.012, iKsl = 0.996 ±0.013, = 0.0413 ± 0.0015, (19) 

\V;t,Vtd\ = 0.0083 ± 0.0016 \VtbV*\ = -0.047 ± 0.008, 

and the CKM phase 7 = 60° ± 14°, sin(2/5) = 0.736 ± 0.049. 

C. Masses and lifetime 

There are two types of quark mass in our analysis. One type is the pole mass which appears 
in the loop integration. Here we fix them as 

ruu = rud = rus = 0, rric = 1.47 GeV, rub = 4.8 GeV. (20) 
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The other type quark mass appears in the hadronic matrix elements and the chirally enhanced 
factor = =^ through the equations of motion. They are renormahzation scale dependent. 
We shall use the 2004 Particle Data Group data |25 for discussion: 

m„(2Gey) = 0.0015 ~ 0.004 GeV, md{2GeV) = 0.004 ~ 0.008 GeV, 
m,(2Ge1/) = 0.08 ~ 0.13 Gel/, mb(mfe) = 4.1 ~ 4.4 GeV, (21) 



and then employ the formulae in Ref. pi] 



m(/i) = m(/io) 



-,(0) 



1 + 



(22) 



2Po 2[3l ) 47r 

to obtain the current quark masses to any scale. The definitions of /5o, /?i can be found 

in [H]. 

To compute the branching ratio, the masses of meson are also taken from 

= 5.279 GeV, mK'± = 0.892 GeV, mK± = 0.494 GeV, 
niB^ = 5.279 GeV, ruK-o = 0.896 GeV, m^o = 0.498 GeV. 

The lifetime of B meson 



= (1.638 ± 0.011) ps, tb^ = (1.532 ± 0.009) ps. 
D. The LCDAs of the meson 

For the LCDAs of the meson, we use the asymptotic form 



(23) 



$p(a;) =6x(l-x), $J'(x) = l, 



(24) 



for the pseudoscalar meson, and 



$iy(x) = <l>l{x) = g'f^^ = 6x(l - X 



(25) 



for the vector meson. 

We adopt the moments of the $f (0 defined in Ref. pTlj ITT] for our numerical evaluation: 



1 'l'f(0 rriB 
4 Ab 



(26) 
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with Xb = (0.46 ±0.11) GeV [23]. The quantity parameterizes our ignorance about the B 
meson distribution amphtudes and thus brings considerable theoretical uncertainty. 
E. The decay constants and form factors 

For the decay constants, we take the latest light-cone QCD sum rule results (LCSR) jT2] in our 
calculations: 

=0.161 Ger, /k = 0.160 Ger, = 0.217 GeV^, = 0.156 Gel^. (27) 

For the form factors involving the B K^*"* transition, we adopt the the values given by ^3] 

0.374 ± 0.034, Af"''>^'^*(0) = 0.292 ± 0.028, 
0.259 ± 0.027, (0) = 0.411 ± 0.033, (28) 

0.331 ±0.041. 

4 Numerical results and Analysis 

First, we will show our estimations in the SM by taking the center value of the input parameters 
and compare with the relevant experimental data. Then, we will consider the RPV effects 
to constrain the relevant RPV couplings from the experimental data. Using the constrained 
parameter spaces, we will give the RPV SUSY predictions for the branching ratios, the direct 
CP asymmetries and the longitudinal polarizations, which have not been measure yet in i? — 
system. 

When considering the RPV effects, we will use the input parameters and the experimental 
data which are varied randomly within la variance. In the SM, the weak phase 7 is well 
constrained, however, with the presence of the RPV, this constraint may be relaxed. We 
would not take 7 within the SM range, but vary it randomly in the range of to vr to obtain 
conservative limits on RPV couplings. We assume that only one sfermion contributes at one 
time with a mass of 100 GeV. As for other values of the sfermion masses, the bounds on the 
couplings in this paper can be easily obtained by scaling them with factor = ( loocey 

For the B Ki*)K^*^ modes, several branching ratios and one direct CP asymmetry have 
been measured by BABAR, Belle and CLEO (211 12^], and their averaged values j2ZI are 

B{B+ K+K") = (1.2 ± 0.3) X 10"^ 
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A^-'''^''\0) = 



^ K^K^) = {OMt^oiD X 10-^ 

B{B+ ^ J^+i?*°) < 5.3 X 10-^ ( 90% CL ), 

B{B+ K*+K*^) < 71 X 10"^ ( 90% CL ), 

^ i^*Oi?*0) < 22 X 10"*^ ( 90% CL ), 

A^jpiB^ ^ K+K^) = 0.15 ±0.33. (29) 

The numerical results in the SM are presented in Table I, which shows the results for the 
CP averaged branching ratios (B), the direct CP asymmetries (Aqp) and the longitudinal 
polarization fractions (/l)- 

Table I: The SM predictions for B (in unit of lO'^), A'^^p and 
in i? ^ K^*^K^*^ decays in the framework of NF and QCDF. 



Decays 


B 


Adir 


Jl 


NF 


QCDF 


NF 


QCDF 


NF 


QCDF 


B+ ^ J^+i?o 


0.61 


0.89 


0.00 


-0.13 








0.57 


0.89 


0.00 


-0.13 






B+ K*+K° 


0.06 


0.10 


0.00 


-0.19 






B+ K+K*° 


0.15 


0.18 


0.00 


-0.08 






B^ K*^K^ 


0.05 


0.10 


0.00 


-0.18 






B^ K^K*^ 


0.14 


0.16 


0.00 


-0.10 






B+ K*+K*° 


0.20 


0.22 


0.00 


-0.22 


0.91 


0.90 


B^a K*^K*° 


0.19 


0.20 


0.00 


-0.22 


0.91 


0.90 



From Table I, we can see that the branching ratios for them are expected to be quite small, of 
order 10~ , since are the pure b ^ d penguin dominated decays. The subleading 

diagrams may lead to the significant CP violations in the most B decays. As 

B^ — * K^K^ decays involved only non-factorizable annihilation contributions, their branching 
ratios are much smaller than those of i? ^ K^K^, K^K^ decays, we would not study the 
B^ K^K^ modes in this paper. It should be noted that the amplitude for B^ K^K*^ is 
not simply related to that for B^ K^K*^ since the spectator quark is part of the in the 
latter decay, while in the former in the K*^. 
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Although recent experimental results in B ^ seem to be roughly consistent with 

the SM predictions, there are still windows for NP in these processes. We now turn to the RPV 
effects in B ^ }{(*)}{(*) decays. There are five RPV coupling constants contributing to the 
eight B Ki*)K^*^ decay modes. We use B, A^p and the experimental constraints shown in 
Eq. ()29|) to constrain the relevant RPV parameters. As known, data on low energy processes 
can be used to impose rather strictly constraints on many of these couplings. In FigEl we 
present the bounds on the RPV couplings. The random variation of the parameters subjecting 
to the constraints as discussed above leads to the scatter plots displayed in Fig0 




Figure 5: The allowed parameter spaces for the relevant RPV couplings constrained by -B — > 
K^*^K^*\ and (pRpv denotes the RPV weak phase. 

From Figl^l we find that every RPV weak phase has two possible bands, one band is for 
positive value of RPV weak phase, and another for negative one. We also find the magnitudes 
of the relevant RPV couplings have been up limited. The upper limits are summarized in Table 
II. For comparison, the existing bounds on these quadric coupling products IHH] are also listed. 
Our bounds on |A'j]^3A^22l5 |A'ii2Ai32l and |A^23Ai2il stronger than the existing ones. 

Using the constrained parameter spaces shown in FigEl one can predict the RPV effects 
on the other quantities which have not been measured yet in i? ^ /^(*)^(*) decays. With the 
expressions for B, A'^p and J'l at hand, we perform a scan on the input parameters and the new 
constrained RPV coupling spaces. Then the allowed ranges for B, Acp and are obtained 
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Table II: Bounds for the relevant RPV couplings hy B ^ /^(*)/^(*) decays for 
100 GeV sfermions and previous bounds are listed for comparison. 



Couplings 


Bounds [Process] 


Previous bounds [Process] 


1 \" \"* 1 


< 2.9 X 10-3 [B - 


^ _/^(*)_^(*)] 


<5.xl0-3 [B^KK] r^t 
<6.xl0-5 [B0^4>-7TO,4,<j)] |a| 


'^il3'^i22 


< 2.2 X 10-3 [B - 


^ _/^(*)_^(*)] 


< 2.9 X 10-^[B KK] m 


'^i22-^'i31 


< 1.7 X 10-3 [B - 




< 1. X 10-^ [KK]^ 


'^'il2-^i32 


< 3.0 X 10"^ [B - 




< 4. X 10"^ [5° ^ 07r°] m 


■^'123-^^1 


< 3.0 X 10-^ [B - 




< 4. X 10"^ [5° ^ 07r°] m 



Table III: The theoretical predictions for B (in unit of 10 A'^p and /l base on the RPV 
SUSY model, which are obtained by the allowed regions of the different RPV couplings. 





\ // \ //* 

'^j23'^il2 


'^il3'^i22 


^^422^^31 


'^il2'^i32 


'^i23'^i21 


BiB+ K*+K°) 


[0.0052, 7.8] 


[0.013,5.5] 


[0.0059,6.4] 


[0.056, 1.4] 


[0.064, 1.3] 


B{B+ K+K*^) 


[0.071,5.3] 


[0.056,5.3] 


[0.0096, 5.3] 








[0.0060, 7.5] 


[0.011,5.1] 


[0.0053,6.1] 


[0.049, 1.5] 


[0.054, 1.2] 


B{B^^ K^K*^) 


[0.069,5.0] 


[0.050,5.1] 


[0.0093,5.0] 






B{B+ K*+K*^) 


[0.087, 19] 


[0.041,23] 


[0.029, 16] 






B{B^ K*^K*^) 


[0.080, 17] 


[0.039, 22] 


[0.027, 15] 






A'^^piB^, - K^K^) 


[-0.75,0.57] 


[-0.19,0.44] 


[-0.18,0.47] 


[-0.18,0.47] 


[-0.18,0.50] 


■^CpiB+ K*+K°) 


[-0.19,0.17] 


[-0.32,0.17] 


[-0.42,0.47] 


[-0.99,0.99] 


[-0.98,0.76] 


A'^UK K+K*^) 


[-0.63,0.63] 


[-0.38,0.47] 


[-0.65,0.38] 








[-0.28,0.19] 


[-0.33,0.17] 


[-0.28,0.80] 


[-0.99,0.99] 


[-0.99,0.73] 




[-0.76,0.62] 


[-0.38,0.48] 


[-0.39,0.40] 






Af^'piB^ K*+K*^) 


[-0.63,0.30] 


[-0.26,0.25] 


[-0.77,0.32] 








[-0.46,0.38] 


[-0.26,0.25] 


[-0.77,0.32] 






fL{B+ ^ K*+K*'') 


[0.72,0.97] 


[0.59,0.95] 


[0.74,0.93] 






fUBj ^ K*0K*^) 


[0.72,0.97] 


[0.59,0.95] 


[0.74,0.93] 
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with five different RPV couplings, which satisfy all present experimental constraints shown in 

We obtain that the RPV effects could alter the predicted B and Aq<p significantly from their 
SM values. For decay modes, which have not been measured yet, their branching ratios can be 
changed one or two order(s) of magnitude comparing with the SM expectations, 

9. X 10"^ < B{B K+K*\ K^K*^) < 5. x 10~^ 
5. X 10"^ < B{B K*+K^, K*^K°) < 8. x 10-^ 
3. X 10"^ < B{B K*+K*^, x*°x*o-^ < 2. x 10-^ (30) 

especially, the upper limit of B{B — K*^K*^) < 2. x 10^^ which we have obtained is smaller 
than the experimental upper limit < 7. x 10^^. For Aqp, the RPV predictions on two decays 
B K*+K*^, g^rg 

A'^SUB IC+K*^) < 0.32, ^cp(^ ^ < 0.38, (31) 

and there are quite loose constraints on the direct CP asymmetries of the other five decays B — >• 
K^K^, K*+K°, K+K*^, K*°K^, K^K*^. But the RPV effects on the hiB K*+K*^, x*^x*o^ 
are found to be very small, fiiB j{*o ^ j^*o j^*o^ found to lie between 0.7 and 1, 

and these intervals are mainly due to the parameter uncertainties not the RPV effects. So 
we might come to the conclusion, the RPV SUSY predictions show that the decays B — >■ 
K*'^K*^, K*^K*^ are dominated by the longitudinal polarization, and there are not abnormal 
large transverse polarizations in Bu^ K*K* decays. The detailed numerical ranges which 
obtained by different RPV couplings are summarized in Table III. 

In FigslUHTUl we present correlations between the physical observable B, A^p, Jl and the 
parameter spaces of different RPV couplings by the three-dimensional scatter plots. The more 
information are displayed in FigslHHTIH we can see the change trends of the physical observable 
quantities with the modulus and weak phase of RPV couplings. We take the first plot in 
Figini as an example, this plot shows that B{B K*^K^) change trend with RPV coupling 
A^23'^ii2- We also give projections on three vertical planes, the \)^'l2'iKi2\'4'ijip plane display the 
allowed regions of y^i2z^ii2 which satisfy experimental data in Eq. ()29|) (the same as the first plot 
in FigEI). It's shown that B{B K*+K°) is increasing with |A.23A.;*2| on the B{B K*+K^) 
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-\K2zKi2\ plane. From the B{B K*+K^)-(j)i^ plane, we can see that B{B K*+K^) is 
increasing with 10^^!. Further refined measurements of B{B K*'^K'^) can further restrict 
the constrained space of A^23Afi*2! whereas with more narrow space of A"23A"*2 more accurate 
B{B K*+K°) can be predicted. 

The following salient features in Figs IMTOl are summarized as following. 

• Figl displays the effects of RPV coupling A^'23A-i2 on B, A'^i'p and in 5 K'^*'>K^*\ 
The constrained \^i23Ki2\~^^p plane shows the allowed range of \i23^ii2 as in the first 
plot of FigEl The six B{B K*+ K'^ , K+ K*^ , K*'^ K'^ , K*'^ , K*+ , K*'^ K*'^) have 
the similar change trends with |A'/23Afj*2l and \(f)^^ |, and they are increasing with |A'/23Afj*2l 
and |. IAqp^B K'^K^)\ are increasing with 10^^ |, but |A'/23A'/j*2| has small effect on 
A'^'piB K^K^). The two \A'}i'p{B K+ K*^ , K*^)\ tend to zero with increasing 
l'^'/23^'/i2l and |0^J. The other four \A'^SpiB K*+ , K*^ K° , K*+ K*^ , K*^ K*^)\ tend 
to zero with increasing \(f)^^ |, and they could have smaller ranges with |A'/23Afi2l- The RPV 
effects on the ^ k*+ K*'^ , K*^ K*'^) are very small, and ^ k*+K*^, K*^K*'^) 
are found to lie between 0.72 and 0.97. 

• The effects of Xii3X'i*22 on B, Acp and fi are exhibited in Fig0 The constrained |A^]^3A'j22|- 
(pj^p plane is the same as the second plot in Fig El The effects of A'ji3A'j22 on B, A^p and 
fi are similar to Af23Afi2 shown in FigEl 

• In FiglHl we plot B, A^p and fi as functions of A^22Ai3i- The constrained |A^22A'i3i|- 
(f)^p plane is the same as the third plot of FigEl The six branching ratios are in- 
creasing with |A^22'^'t3il and decreasing with \Aqp{B — > K^K'^)\ is unaffected 
by |A'j22A^3i|, but the other six direct CP asymmetries could have smaller ranges with 
|Ai22Ai3il- lAcpiK*"^ K° , K*'^ K'^)\ tends to zero with decreasing l^^^l, however, has 
small effect on Af^^B ^ K°K°, K+K*^, K'^K*^, K*+K*^, K*°K*°). The A^22A-3i effects 
on the hiB k*+ K*'^ , K*'^ K*°) are small. 

• RPV coupling A'^i2A-*32 contributes to the decays B K+K°, K^K^, K*+K^, K*^K^, and 
the effects are shown in FiglHl The constrained |A^]^2Ai32|-09p plane is the same as the 
fourth plot in Fig0 We can see that B{B K*+K^,K*^K^) are rising with |A'ii2A-32|, 
and unaffected by . Acp{B K^K^) is steady against |A-i2A-32|, and \Acp{B — > 



K*+K^,K*'^K°)\ could have smaller ranges with {X'a^^^^l- ^cp(^ ^ K^K°,K*+K^, 
K*^K^) are becoming large with increasing of 10^^ |. 

• ^i23^i2i also only contributes to the decays B K+ K'^ , K° , K*+ , K*^ , and its 
effects are shown in FigEl The constrained |A^23-^'i2il"09p plane is the same as the last 
plot in Figini B{B K*~^ , K*^ K^) are increasing with |A-23A-2i|, and unaffected by 
-^cpiB K^K^) is steady against lA'^ssA^il. and \Af.'p{B K*+ , K*^ K^)\ 
could be varied in small ranges with |A^23A-2il- ^cp(^ K°K°) is decreasing with 
109^ I, but Afi'piB K*+K^, K*^K^) are increasing with 10^^ |. 

The predictions of B and are quite uncertain in the RPV SUSY, since we just have few 
experimental measurements and many theoretical uncertainties. One must wait for the error 
bars to come down and more channels measured. With the operation of B factory experiments, 
large amounts of experimental data on hadronic B meson decays are being collected, and 
measurements of previously known observable will become more precise. From the comparison 
of our predictions in Figs IMTOl with the near future experiments, one will obtain more stringent 
bounds on the product combinations of RPV couplings. On the other hand, the RPV SUSY 
predictions of other decays will become more precise by the more stringent bounds on the RPV 
couplings. 

5 Conclusions 

In conclusions, the pure penguin B decays are very important for understanding 

the dynamics of nonleptonic two-body B decays and testing the SM. We have studied the 
B K^*)k'^*) decays with the QCDF approach in the RPV SUSY model. We have obtained 
fairly constrained parameter spaces of the RPV couplings from the present experimental data 
of B ^ I&)K'^*) decays, and some of these constraints are stronger than the existing ones. 
Furthermore, using the constrained parameter spaces, we have shown the RPV SUSY expec- 
tations for the other quantities in i? ^ /('(*)^(*) decays which have not been measured yet. 
We have found that the RPV effects could significantly alter B and from their SM values, 
but fL{B K*^K*^, ]^*0]^*o^ g^j,g T^Q^ significantly affected by the RPV effects and the decays 
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B — > K*~^K*^, ^*o^*o g^jjj dominated by the longitudinal polarization. We also have pre- 
sented correlations between the physical observable B, Af}p, /l and the constrained parameter 
spaces of RPV couplings in FigslHlfTIH which could be tested in the near future. 
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Appendix 



A. Correction functions for B M1M2 decay at order 

In this appendix, we present the explicit form for the correction functions appearing in the 
parameters and b^. It's noted that in i? ^ PV decays, $a/(m) "^if ('^) if M a vector 
meson. 

A.l The correction functions in B ^ PP,PV decays 

• One-loop vertex correction function is 

TTih /I — 2,u \ 

= 121n — - 18 + 3 / dui Inu-iTr] $M,(n). (32) 

u Jo \ 1 — u J 



• The hard spectator interactions are given by 

'''''''' ~ mlFr'Hmli.) Jo T*^ ^^Vo T^^^^^^Vo T 
if M2 is a pseudoscalar meson, and 



U Jo V 



,(33) 



if M2 is a vector meson. 

Considering the off-shellness of the gluon in hard scattering kernel, it is natural to associate 



a scale /i^ ~ y ^QCD'^b , rather than /i ~ m^. For the logarithmically divergent integral, we will 
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parameterize it as in [T7]: Xh = du/u = —ln{!s.QCD/^b) + Qh^^'^" ^b/ ^qcd with {qhi4'h) 
related to the contributions from hard spectator scattering. In the numerical analysis, we take 
^QCD = 0.5GeV, {gh, (pn) = (0, 0) as our default values. The same as in i? — VV decay. 
• The penguin contributions at the twist-2 are described by the functions 



pp 



M2,2 



CiGm,{sp) + Gs GmM + Gm.(1) + (G, + G,) (uf - 2)GmM 



2nf 



PfifS" = {Ci + NcG,)Gm.As,)- G^.^f J^du^ 



'0 1 — u 
1 , 3$m2(m) 



u 



(35) 



where rij = 5 is the number of quark flavors, and s„ = 0, Sc = (mc/mfe)^ are mass ratios involved 
in the evaluation of the penguin diagrams. The function ^^/^(s) is defined as 



Gm2(s) = - + - In 1-4 du dx XX In (s — xxM — ?e)$j\/2('w)- (36) 

3 3 /i Jo Jq 

• The twist-3 terms from the penguin diagrams are given by 



pp 



pP,EW 



G,Gm2{sp) + Gs Gm2(0) + Gm,(1) + (C4 + C6) 



inj-2)GMM 



GKui^cj + ^4/2(1 



{G^ + NcG2)Gm2{sp)- -iG^ 



(37) 



with 



1 /■! 



Gm2(s) = ^ + ^ln— + 4 \' du r dx xxlnis - xxu - ie)^^^Hu), (38) 
3 3 fi Jo Jo ^ 

if M2 is a pseudoscalar meson, and we omit the twist-3 terms from the penguin diagrams when 
M2 is a vector meson. 

• The weak annihilation contributions are given by 



18 I -4+ -) + 2rfrfXl 



Al{M,, M2) ^ 67ra,(rf ^ - rf ^) ( - + — ) , 



vr 



A{(Mi, M2) ^ 67ra,(rf + rf ^) (2X1 - Xa 
A((Mi,M2) = 0, 4(Mi,M2) = 0. 



(39) 
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when both final state mesons are pseudoscalar, whereas 

Ai(Mi, M2) ^ -v4^(Mi, M2) ^ ISvra, f - 4 + 



TT 



2' 



A'3(Mi, M2) ^ 67ra,r^^^^ (^X^ - 2Xa + yj , 

4(Mi, M2) ^ -evra.rf (2X1 " Xa) , 

M2) = 0, M2) = 0. (40) 

when Ml is a vector meson and M2 is a pseudoscalar. For the opposite case of a pseudoscalar 
Ml and a vector M2, one exchanges r^^ ^ r^^ in the previous equations and changes the sign 
of 4. 

Here the superscripts i and / refer to gluon emission from the initial and final state quarks, 
respectively. The subscript k of A^jl^ refers to one of the three possible Dirac structures Fi ® F2, 
namely k = I loi {V - A)®{V - A), k = 2 for {V -A)®{V + A), and A; = 3 for {-2){S-P)® 
(S + P). Xa = Jq dn/n is a logarithmically divergent integral, and will be phenomenologically 
parameterized in the calculation as Xh- As for the hard spectator terms, we will evaluate the 



various quantities in Eqs. (jHIIjl and at the scale fih = y-^QCDfnt,. 
A. 2 B ^VV decays 

In the rest frame of B system, since the B meson has spin zero, two vectors have the same 
helicity therefore three polarization states are possible, one longitudinal (L) and two transverse, 
corresponding to helicities A = and A = ± ( here Ai = A2 = A). We assume the Mi(M2) meson 
fiying in the minus(plus) z-direction carrying the momentum pi(p2), Using the sign convention 

g0123 ^ Yiaye 



A 



M1M2 



ho, 



(41) 



where for A = and h± for A = ± 

nih , „ , \ 1 ~ 2m 



• V^A (±1) contain the contributions from the vertex corrections, and given by 



1/° (a) = 121n— - 18+ / du^fHu) [3 Inw-SiTr , (42) 

H JO " V 1 — M 



l/4(a) = 121n^-18+ (" du U^'Hu) fsl^ In n - 3zvr^ 

jji Jo \ A j \ 1 — u J 
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• The hard spectator scattering contributions, exphcit calculations for H^j^j^j^{a) yield 



Nc ho Jo ^ JO V Jo u 



47r^ 2ifBfMjM2mM, 
Nc mBh± 
1 



Jo t Jo 



dv- 



y2 



l(a)M2 I \ 

X / du\^^''^u)-'^^^ 







X / dt^dn 
Jo 



(v)M: 



4 

/(a)Mi 



e Jo 

4:71^ 2ifB /m, /m2 rriM, ^^^liO 
Nc m%h± Jo ^ 



V) \ I (v)M2/ \ , ^' 



l{a)M2 



{u) \ U + V 



(43) 



with V = 1 — V, when the asymptotical form for the vector meson LCDAs adopted, there will 
be infrared divergences in i?MiA4'2- OEHl, we introduce a cutoff of order kqcolmh and 

take IS.QCD = 0.5 GeV as our default value. 

• The contributions of the QCD penguin-type diagrams can be described by the functions 



^M2,2 



pX,p,EW 
^M2,2 



q—u 



+\Co[e,Glj^{s,) + e^G\,^{s,)\ + ^(Cg + G^o) E e<?' 



2 '-'77 '-^g' 



(44) 



G^2(«) = ^ + ^ln^ + 4 /'dn^f (u)g(u,.). 



2 2, mfc /-i 
- + -ln^ + 2 / dM 

O O /i JO 



l{a)M2 



U 



with the function g(u, s) defined as 

g(u, s) 



dx xa; In is — xxu — ie) 



{u,s), 



(45) 



(46) 



We omit the twist-3 terms from the penguin diagrams for B VV decays. 

• We have also taken into account the contributions of the dipole operator Og^, which are 
described by the functions 

1 2^1'' (u) 







- du 




Jo 




/■I du 




Jo u 



1-u 



-ugy [u) T 



UK ' Uu) 



+ dv[^f'{v) 



{y)M2 



(a)M2 



U 



(47) 
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here we consider the higher-twist effects /c'* = uEnt + kj_ + -^n'^ in the projector of the vector 
meson. The = in Eq. (fTTjl j2Sl IHO] if considering the Wandzura- Wilczek-type relations 

I2n!. 

We have not onsidered the annihilation contributions in B VV decays. 
A. 3 The contributions of new operators in RPV SUSY 



Compared with the operators in the 'H^^, there are new operators {q2(l3)v±A{bqi)v+A in the 
• For B PP, PV decays, since 



{P\ gi7/.(l - 75)^2! 0) = -{P\ 517^(1 + 75)92! 0) = -{P\ gi7/.75g2| 0), 
(P| g7^(l - 75)6 I B) = {P\ g7^(l + 75)6 | B) = {P\ q^,b \ B), 
{V\ qil,{l-lM 0) = {V\ qil,{l+lM 0) = {V\ gi7^g2| 0), 
{V\ q^,{l - 75)6 I B) = -{V\ 57^(1 + 75)& I B) = -{V\ q^,^,b \ B), (48) 

the RPV contribution to the decay amplitude will modify the SM amplitude by an overall 
relation. 

• For B — > VV, we will use the prime on the quantities stands for the {q2q3)v±A(bqi)v+A 
current contribution. In the NF approach, the factorizable amplitude can be expressed as 

^M.M, = (M2|fe7M(l-«75)g3)|0)(Mi|(67^(l + 75)gi)|5). (49) 
Taking the Mi(M2) meson flying in the minus(plus) z-direction and using the sign convention 

g0123 ^ _1 \]_Q^^Q 



A' 



(ml - ml,^ - minims + m^JAf^^^^ {ml, J - ^^^^A^^'^^ {ml, J 



(50) 



The vertex corrections V^}^{a) and the hard spectator scattering corrections H'^^j^,j^{a) as 
follows: 

v't^iAa) = -121n— + 18-6(l + a)- /'du$f^(u) fs^^^lnu-SiTrV 

Vl±(a) = -12 In ^ + 18 - 6(1 + a) - du Ul>''\u) ± "g^"^^""^ (") (^3!^ \nu- 3..) , 
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Nc K Jo i Jo V Jq 

""-^-^^"^ = ^'^'Vo '^^io '"^^ 



A ( (v)M2( \ ag[^ {u)\ Att^ 2tfBfMifM2mMimM2 f ^c^iiO 
dulgy {u) + +— / — ^ 

X / d.d. g^l^^'^iv) t'-^ ^ ) g'l'''^{u) ±^ ^ 1 



B . The amplitudes in the SM 



K°K*° 



04 - -aio + r 
2 ^ 



(ae - -as) 



V2 



tb^td 



1 K°, 1 N 

a4--a.io + r^ [ag - -as) 



} 



-Vtiytd 



V2\ 



tb^td 



b3{K", K") + b4{K", K") + b4(K", K") 



-^bf"iK'^,K°) - ^bf^(K°,K°) - ^bl"' {K° , K'^) 



a4 - -aio - r. 



(ae - -as) 



^fBfKfK' {v:,V^db2iK*+,K0)-V,lVta [b3{K*+ , K°) + b^-^ {K*+ , RO)]} 



tb^td 



14, - -aio 



} 



V2 

^fsfKfK' {v:,V„ab2{K+,K''') - v:,Vu [63(i^+,i?*°) + f)r(/f+, A-*")] } : 
V2 



{-Vtl 



tb^td 



a4 - -aio - ("6 - -as) 



^fBfKfK' {- 



bi{K'°, K°) + b4{K*°,K°) + b4{K°, K'°) 
1 . 



.0 ^?0^ 



2 

G F 

7f 



} 



^tb^td 



04 - -aio 



b-i{K°,K*'^) + b4{K°,K*°) + b4{K*°, K°) 



■hl'^{K°,K*°) - if)|"'(A'°,_ft:*0) - ^bl'^ {K*° , K°) I 



tb^td 



a4 - -aio 



a4 — 2'^10 



Here we have not considered the annihilation contributions in i? — VV decays. 



C . The ampHtudes for RPV 

-■■ ( Ml Ml* 

A'^ p {B+ ^ K+ K°) = 



16m? 



\ll \//* / \/ \/* \f \ 



8m? 



r+ro 
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^^fp(BO ^ K°K°) 
A^p(B+ K*+K°) 

A'^p(B+ — K+K*°) 
^^p(bO ^ K*°K°) 

A'^p{B+ K*+K*°) 



8m? 



8m? 



8m? 

"Li 



8m? 



-8//3o if^ 

X 



8m? 

''Li 

X 



8m? 



16m?, 



-4/00 ^ 



8m? 



^il2'^»32 

8m? 



^i23^i21 

8m? 



X 



16m? 
^i23 il2 . 



K+K*0 



16»n? 



-4/00 ^ 



'^il2^»32 
8m? 



^i23^i21 
8m? 

•^Li 



^i23^il2 „-4/0o P - _i_ 



16m? 



_i23_il2_ 
16m?. 



-4/00 p 



8m? 

\ / \/* 

^il3 i22 

8m? 

I^Li 

-8/00^^0 
X 



8m? 

\/ \ /* 
^113 t22 

8m? 



^i22^i31 
8m? 

I^Li 



Sm'i 

•^Li 

\f \ 

i22^i31 

8m? 



-8/00 2, 



K+K'O 



^-8/0O(_i^.„^„) 



-8/00 r' /l' 



A'fp{B+ K*°K*°) 



8m? 



^»23^il2 -4/00 p/ 



16m? 



_i22_m_ 

8m? 



^il3^i22 
8m? 



-8/00 i 



A' 



K'OK'O 



In the ^^J', F^j^j^i^ and Lm^^.j^ are defined as 
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LmiM2 

PP, PV decays, and 



^ ]_ as Cf 

Nc 47r A^c 

A^c I 47r Nc 



VkPi + HmiM2 

12 + Va/2 + HmiM2 



F' 



Tl 

^MiM2 



1 - 


1 








AtiNc 


1 




Us Cf 




{- 


AtiNc 


1 




as Cf 




{- 


A-K Nc 



ylX 



i2 + v^Xi^) + H''m 



i2 + v^4( 



-l)+^^M.A/2(-l) 



B — > VV decays. 
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Figure 6: The effects of RPV coupling AfsgAffg in B ^ k^*)k^*) decays. 
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Figure 7: The effects of RPV coupling A'^^gAfsa in 5 ^ k'^*)k'^*) decays. 
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Figure 8: The effects of RPV coupling ^22^^,1 in -B ^ k'^*)k'^*) decays. 
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Figure 9: The effects of RPV coupling A'^igAfga mB ^ k^*)k'^*) decays. 
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1 Introduction 



The study of exclusive hadronic S-meson decays can provide not only an interesting avenue to 
understand the CP violation and flavor mixing of the quark sector in the standard model (SM), 
but also powerful means to probe different new physics (NP) scenarios beyond the SM. Recent 
experimental measurements have shown that some B decays to two light mesons deviated from 
the SM expectations, for example, the tttt, ttK puzzle P and the polarization puzzle in B VV 
decays [2]. Although these measurements represent quite a challenge for theory, the SM is in no 
way ruled out yet since there are many theoretical uncertainties in low energy QCD. However, 
it will be under considerable strain if the experimental data persist for a long time. 

Among those NP models that survived electroweak (EW) data, one of the most respectable 
options is the R-parity violating (RPV) supersymmetry (SUSY). The possible appearance of 
the RPV couplings [Hj, which will violate the lepton and baryon number conservation, has 
gained full attention in searching for SUSY [HE]. The effect of the RPV SUSY on B decays 
have been extensively investigated previously in the literatures IZj ? and it has been proposed 
as a possible resolution to the polarization puzzle and the vrvr, nK puzzle |H]. The pure penguin 
B K^*^K^*^ decays are closely related with the puzzles which are inconsistent with the SM 
predictions, and therefore are very important for understanding the dynamics of nonleptonic 
two-body B decays, which have been studied in Refs. jHj. If the RPV SUSY is the right model 
to resolve these puzzles, the same type of NP will affect B k(*)k^*^ decays. In this work, we 
shall study the RPV SUSY effects in the B ^ K^*)Ki*) decays by using the QCD factorization 
(QCDF) approach jTUj for hadronic dynamics. The B }{(*)}{(*) decays are all induced at 
the quark level by 6 ^ dss process, and they involve the same set of RPV coupling constants. 
Using the latest experimental data and the theoretical parameters, we obtain the new upper 
limits on the relevant RPV couplings. Then we use the constrained regions of parameters to 
examine the RPV effects on observations in the B — > K^*^K^*^ decays which have not been 
measured yet. 

The paper is arranged as follows. In Sec. 2, we calculate the CP averaged branching ratios, 
the direct CP asymmetries oi B ^ and the polarization fractions in i? K*K* 

decays, taking account of the RPV effects with the QCDF approach. In Sec. 3, we tabulate 
the theoretical inputs in our numerical analysis. Section 4 deals with the numerical results. 
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We display the constrained parameter spaces which satisfy all the experimental data, and then 
we use the constrained parameter spaces to predict the RPV effects on the other observable 
quantities, which have not been measured yet in B system. Section 5 contains our 

summary and conclusion. 



2 The theoretical frame for B K^*^ K^*^ decays 
2.1 The decay amplitudes in the SM 

In the SM, the low energy effective Hamiltonian for the AB = 1 transition at the scale n is 
given by [TT] 



r 10 ^ 

'^e// ~ ~~/^ X/ ^P] ^^Ql + ^2^2 + \CiQi + Cj^Qj^ + CsqQsq >+h.C. 
V ^ p=u,c [ 1=3 ) 



(1) 



here \p = VphV*^ for 6 — > g transition (p G {m, c}, q G {rf, s}) and the detailed definition of the 
operator base can be found in [TT] . 

Using the weak effective Hamiltonian given by Eq.Q, we can now write the decay ampli- 
tudes for the general two-body hadronic B — > M1M2 decays as 

A^^{B ^ M1M2) = {MiM2\H^fj\B) 

= %EE^C'.(/i)(MiM2|Q,(/x)|E). (2) 

The essential theoretical difficulty for obtaining the decay amplitude arises from the evaluation 
of hadronic matrix elements {MiM2\Qi{^)\B) . There are at least three approaches with different 
considerations to tackle the said difficulty: the naive factorization (NF) ^21^], the perturbative 
QCD [m, and the QCDF The QCDF developed by Beneke, Buchalla, Neubert and 

Sachrajda is a powerful framework for studying charmless B decays. We will employ the 
QCDF approach in this paper. 

The QCDF ^01 allows us to compute the nonfactorizable corrections to the hadronic matrix 
elements {MiM2\0i\B) in the heavy quark limit. The decay amplitude has the form 



(3) 
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Figure 1: The next to leading order nonfactorizable contributions to the coefficients a^. 




Figure 2: The weak annihilation contributions to the coefficients 6f . 



where the effective parameters af including nonfactorizable corrections at order of a^. They are 
calculated from the vertex corrections, the hard spectator scattering, and the QCD penguin 
contributions, which are shown in FigH The parameters 6f are calculated from the weak 
annihilation contributions as shown in Fig|21 

Under the naive factorization (NF) approach, the factorized matrix element is given by 



A 



M1M2 



^ (M2|fe7/.(1 - l,)q-m{M,\{h^{l - 75)gi)|5). 



(4) 



In term of the decay constant and form factors [T3], AmiM2 ^-^e expressed as 



A 



M1M2 



2 T7B-*Mi, 







m 



M2/' 



-ifM^mM^ 



2 AB~*Ah, 







m 



if Ml = P, M2 = P, 
if Ml = P, M2 = V, 
if Ml = 1/, M2 = P, 



{el ■ sDimn + m^JAf ^""Hmi.J - (el ■ pB)ie*, ■ pb) 



2A. 



+^e^lua/3£2^l PbPi 



2^ 



if Ml = V, M2 = V, 



(5) 



where P(V) denote a pseudoscalar (vector) meson, Pb{^b) is the four-momentum(mass) of the 
B meson, 771^. is the masses of the Mj mesons, and e* is the polarization vector of the vector 
mesons M,-. 



Following Beneke and Neubert (TH], coefficients af can be split into two parts: af = of j+af 
The ffist part contains the NF contribution and the sum of nonfactorizable vertex and penguin 
corrections, while the second one arises from the hard spectator scattering. The coefficients 
read [15] 



as,/ 



a 



4J 
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Nc 
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Nc 



1 + 
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1 + 
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0^2,11 



C2 


Cpas 
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Ci 




Nc 


An 




Cpas 



+ 
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An 
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c 



6 ■ 



Cpas 
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Nm, + 
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An 



-12 - V, 



M2, 



Cs 
■ C*io 



1 -6 ■ 



Cpas 



An 



Nm, + 



An Nc 



r,P,EW 

9n Nc ' 



0-A,II 
0-5,11 
06,11 
07,11 



Nc An 

C3 Cpa 



Nc An 

Nc An 
0, 



HmiM2i 
HmiM2) 

HmiM21 

—HM1M2] 



J- H : 



An 
An 



Nc An 
as,// = 0, 



'-H 



A/1M2J, 



V, 



M2 



+ 



Up 



-yP,EW 
A/2,2 



97r N. 



09,11 - 
aio,// 



c 



Nc An 
Cq CpOi. 
' Nc 47r 



HmiM2, 

-HmiM2j (6) 



where = as(//), Cp = {N^ — l)/{2Nc), Nc = 3 is the number of colors, and Nm2 = 1(0) for 
M2 is a pseudoscalar (vector) meson. The quantities Vm2, -f^A/iA/2) P^h 25 ^Ah 3' PhifY PmT 
consist of convolutions of hard-scattering kernels with meson distribution amplitudes. Specifi- 
cally, the terms Vm2 come from the vertex corrections in Figma)-md), -Pm2,2 -^A/2,3 (P^h^ 
and Pm^Y) arise from QCD (EW) penguin contractions and the contributions from the dipole 
operators as depicted by FigHl^e) andlH^f). HmiM2 is due to the hard spectator scattering as 
Fig^g) and^h). For the penguin terms, the subscript 2 and 3 indicate the twist 2 and 3 
distribution amplitudes of light mesons, respectively. Explicit forms for these quantities are 

relegated to Appendix A. 

We use the convention that Mi contains an antiquark from the weak vertex, for non-singlet 
annihilation M2 then contains a quark from the weak vertex. The parameters 6f = 6f (Mi, M2) 
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in Eq.Q correspond to the weak annihilation contributions and are given as ^7] 



&i(Afi, M2) = -^Ci^KMi, M2), ?>2(Mi, A/2) = -^C2Ai(Mi, M2), 

&f (Ml, M2) = [Cs^KMi, M2) + C5 (^^(Mi, M2) + 4(Mi, M2)) + NcC^Ai{Ah,M2) 
6?(Afi, M2) = [CMNh,M2) + C6^^(Mi, A/2)] , 

6f'^'"(Afi, Af2) = [CgAUAfi, A/2) + C7 (a^(A/i, Afs) + A^(Afi, Afs)) + iVcC84(A^i- ^2) 
65'^'"(Mi, A/2) = -ll^ [CioAl(A/i, A/2) + C8^^(A/i, M2)] , (7) 

the annihilation coefficients (61, 62), 64) and (63''^"', Tf^^^) correspond to the contributions of 
the tree, QCD penguins and EW penguins operators insertions, respectively. The explicit form 
for the building blocks A*/(Mi, M2) can be found in Appendix A. 

With the coefficients in Eq.® and ((Zj), we can obtain the decay amplitudes of the SM 

If 

part) and ^f^^ (the subscript "a" denotes the annihilation part). The SM part amplitudes of 
B K(*)Ki*) decays are given in Appendix B. 



part A^'^^ (the subscript "/" denotes the part without the contribution from the annihilation 



2.2 R-parity violating SUSY effects in the decays 

In the most general superpotential of the minimal supersymmetric Standard Model (MSSM), 
the RPV superpotential is given by 



where L and Q are the SU(2)-doublet lepton and quark superfields and E'^, Lf^ and Z)'^ are the 
singlet superfields, while i, j and k are generation indices and c denotes a charge conjugate field. 

The bilinear RPV superpotential terms fiiLiHu can be rotated away by suitable redefining 
the lepton and Higgs superfields ^Hl- However, the rotation will generate a soft SUSY breaking 
bilinear term which would affect our calculation through penguin level. However, the processes 
discussed in this paper could be induced by tree-level RPV couplings, so that we would neglect 
sub-leading RPV penguin contributions in this study. 

The A and A' couplings in Eq.® break the lepton number, while the A" couplings break 
the baryon number. There are 27 X[j,^ couplings, 9 Xijk and 9 X'-ji^ couplings. X[ij]k are antisym- 
metric with respect to their ffist two indices, and X'-y^ are antisymmetric with j and k. The 

6 



antisymmetry of the baryon number violating couplings A^^^j in the last two indices implies 
that there are no X'l^f. operator generating the b sss and h ddd transitions. 







yi \i* \ X A' \'* 

^ijk ^imn \ / ijk ^imn 

dk dn d„ 

Figure 3: Sleptons exchanging diagrams for nonleptonic B decays. 






A" A"* 



Figure 4: Squarks exchanging diagrams for nonleptonic B decays. 

From Eq.(jH)), we can obtain the following four fermion effective Hamiltonian due to the 
sleptons exchange as shown in FigOl 

id„^YPRdl)s{UkltiPLUj)8, 

{dr,^l^PRdlMdk-i^PLdj)8. (9) 



^2u-2d 


= E 

i 


A' A'* 
2m? . 


-8//3o 




= E 

i 


A' A'* 

^ijm^ikl 

2m? ^ 


-8//3o 



The four fermion effective Hamiltonian due to the squarks exchanging as shown in FigH are 

kh - {ua^PRUj)8{dkj^PRdi)s 
{dkl^PRUj)i{uaf,PRdi)i - {dkYPRUj)8{ua^PRdi)s } , 
{da^PRdj)i{dk-i^,PRdi)i - ida^PRdj)sidk-f^PRdi)s\ ■ (10) 



'^2u~2d 



\" A"* 

E -j^V-'^'''{[iu^rPRuMdkl,PRdi 



A" A"* 



5Z Am.. 



2 ^ 



In Eqs.© and Pl 



1-7; 



1+75 



2 = -^^V = a^m'i,) ^'^'^ /5o = 11 - Inf. The subscript 



for the currents (j^)i,8 represents the current in the color singlet and octet, respectively. The 
coefficients r/~^/^o and ?^~^/^o are due to the running from the sfermion mass scale rrij (100 GeV 
assumed) down to the m?, scale. Since it is always assumed in phenomenology for numerical 
display that only one sfermion contributes at one time, we neglect the mixing between the 
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operators when we use the renormahzation group equation (RGE) to run Tt^jf down to the low 
scale. 

The RPV amplitude for the decays can be written as 

A^{B^MiM2) = (MiM2\7^}f\B) . (11) 

The product RPV couplings can in general be complex and their phases may induce new 
contribution to CP violation, which we write as 

here the RPV coupling constant A G {A, A', A"}, and is the RPV weak phase, which may 
be any value between — vr and vr. 

For simplicity we only consider the vertex corrections and the hard spectator scattering in 
the RPV decay amplitudes. We ignore the RPV penguin contributions, which are expected 
to be small even compared with the SM penguin amplitudes, this follows from the smallness 
of the relevant RPV couplings compared with the SM gauge couplings. The bounds on the 
RPV couplings are insensitive to the inclusion of the RPV penguins [20] • We also neglected the 
annihilation contributions in the RPV amplitudes. The R-parity violating part of the decay 
amplitudes A^^ can be found in Appendix C. 

2.3 The total decay amplitude 

With the QCDF, we can get the total decay amplitude 

A{B M1M2) = A^f^'{B M1M2) + Al^{B ^ M1M2) + A^^ {B ^ MM). (13) 

The expressions for the SM amplitude Aj^ and the RPV amplitude A^^ are presented in 
Appendices B and C, respectively. From the amplitude in Eq. |T!?jl . the branching ratio reads 

}3{B ^ M1M2) = \A{B ^ M1M2) f S, (14) 

where S" = 1/2 if Mi and M2 are identical, and S = 1 otherwise; tb is the B lifetime, |j9c| is the 
center of mass momentum in the center of mass frame of B meson, and given by 

\Pc\ = - i^Mi + "^Af2)^][?Ti| - {niMi - "iA/a)^]- (15) 
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The direct CP asymmetry is defined as 



B{B ^ /) - B{B ^ f) 



B{B^f) + B{B^fy ^^^^ 
In the B VV decay, the longitudinal polarization fraction is defined by 

-n 1/112,1/112,1/112' v-^ ' / 



r 

where Aq{A±) corresponding to the longitudinal(two transverse) polarization amplitude(s) for 
B ^VV decay. 

3 Input Parameters 

A. Wilson coefficients 

We use the next-to-leading Wilson coefficients calculated in the naive dimensional regularization 
(NDR) scheme at mf, scale 

Ci = 1.082, C2 = -0.185, C3 = 0.014, C4 = -0.035, 6*5 = 0.009, 

Cg = -0.041, Cj/ae = -0.002, Cg/tte = 0.054, Cg/a, = -1.292, 

Cio/tte = 0.263, 6*7;^^ = -0.299, = -0.143. (18) 

B. The CKM matrix element 

The magnitude of the CKM elements are taken from [21]: 

\Vud\ = 0.9738 ± 0.0005, |K.| = 0.2200 ± 0.0026, |Kb| = 0.00367 ± 0.00047, 

iKdl = -0.224 ±0.012, iKsl = 0.996 ±0.013, = 0.0413 ± 0.0015, (19) 

\V;t,Vtd\ = 0.0083 ± 0.0016 \VtbV*\ = -0.047 ± 0.008, 

and the CKM phase 7 = 60° ± 14°, sin(2/5) = 0.736 ± 0.049. 

C. Masses and lifetime 

There are two types of quark mass in our analysis. One type is the pole mass which appears 
in the loop integration. Here we fix them as 

ruu = rud = rus = 0, rric = 1.47 GeV, rub = 4.8 GeV. (20) 
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The other type quark mass appears in the hadronic matrix elements and the chirally enhanced 
factor = =^ through the equations of motion. They are renormahzation scale dependent. 
We shall use the 2004 Particle Data Group data |25 for discussion: 

m„(2Gey) = 0.0015 ~ 0.004 GeV, md{2GeV) = 0.004 ~ 0.008 GeV, 
m,(2Ge1/) = 0.08 ~ 0.13 Gel/, mb(mfe) = 4.1 ~ 4.4 GeV, (21) 



and then employ the formulae in Ref. pi] 



m(/i) = m(/io) 



-,(0) 



1 + 



(22) 



2Po 2[3l ) 47r 

to obtain the current quark masses to any scale. The definitions of /5o, /?i can be found 

in [H]. 

To compute the branching ratio, the masses of meson are also taken from 

= 5.279 GeV, mK'± = 0.892 GeV, mK± = 0.494 GeV, 
niB^ = 5.279 GeV, ruK-o = 0.896 GeV, m^o = 0.498 GeV. 

The lifetime of B meson 



= (1.638 ± 0.011) ps, tb^ = (1.532 ± 0.009) ps. 
D. The LCDAs of the meson 

For the LCDAs of the meson, we use the asymptotic form 



(23) 



$p(a;) =6x(l-x), $J'(x) = l, 



(24) 



for the pseudoscalar meson, and 



$iy(x) = <l>l{x) = g'f^^ = 6x(l - X 



(25) 



for the vector meson. 

We adopt the moments of the $f (0 defined in Ref. pTlj ITT] for our numerical evaluation: 



1 'l'f(0 rriB 
4 Ab 



(26) 
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with Xb = (0.46 ±0.11) GeV [23]. The quantity parameterizes our ignorance about the B 
meson distribution amphtudes and thus brings considerable theoretical uncertainty. 
E. The decay constants and form factors 

For the decay constants, we take the latest light-cone QCD sum rule results (LCSR) jT2] in our 
calculations: 

=0.161 Ger, /k = 0.160 Ger, = 0.217 GeV^, = 0.156 Gel^. (27) 

For the form factors involving the B K^*"* transition, we adopt the the values given by ^3] 

0.374 ± 0.034, Af"''>^'^*(0) = 0.292 ± 0.028, 
0.259 ± 0.027, (0) = 0.411 ± 0.033, (28) 

0.331 ±0.041. 

4 Numerical results and Analysis 

First, we will show our estimations in the SM by taking the center value of the input parameters 
and compare with the relevant experimental data. Then, we will consider the RPV effects 
to constrain the relevant RPV couplings from the experimental data. Using the constrained 
parameter spaces, we will give the RPV SUSY predictions for the branching ratios, the direct 
CP asymmetries and the longitudinal polarizations, which have not been measure yet in i? — 
system. 

When considering the RPV effects, we will use the input parameters and the experimental 
data which are varied randomly within la variance. In the SM, the weak phase 7 is well 
constrained, however, with the presence of the RPV, this constraint may be relaxed. We 
would not take 7 within the SM range, but vary it randomly in the range of to vr to obtain 
conservative limits on RPV couplings. We assume that only one sfermion contributes at one 
time with a mass of 100 GeV. As for other values of the sfermion masses, the bounds on the 
couplings in this paper can be easily obtained by scaling them with factor = ( loocey 

For the B Ki*)K^*^ modes, several branching ratios and one direct CP asymmetry have 
been measured by BABAR, Belle and CLEO (211 12^], and their averaged values j2ZI are 

B{B+ K+K") = (1.2 ± 0.3) X 10"^ 
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A^-'''^''\0) = 



^ K^K^) = {OMt^oiD X 10-^ 

B{B+ ^ J^+i?*°) < 5.3 X 10-^ ( 90% CL ), 

B{B+ K*+K*^) < 71 X 10"^ ( 90% CL ), 

^ i^*Oi?*0) < 22 X 10"*^ ( 90% CL ), 

A^jpiB^ ^ K+K^) = 0.15 ±0.33. (29) 

The numerical results in the SM are presented in Table I, which shows the results for the 
CP averaged branching ratios (B), the direct CP asymmetries (Aqp) and the longitudinal 
polarization fractions (/l)- 

Table I: The SM predictions for B (in unit of lO'^), A'^^p and 
in i? ^ K^*^K^*^ decays in the framework of NF and QCDF. 



Decays 


B 


Adir 


Jl 


NF 


QCDF 


NF 


QCDF 


NF 


QCDF 


B+ ^ J^+i?o 


0.61 


0.89 


0.00 


-0.13 








0.57 


0.89 


0.00 


-0.13 






B+ K*+K° 


0.06 


0.10 


0.00 


-0.19 






B+ K+K*° 


0.15 


0.18 


0.00 


-0.08 






B^ K*^K^ 


0.05 


0.10 


0.00 


-0.18 






B^ K^K*^ 


0.14 


0.16 


0.00 


-0.10 






B+ K*+K*° 


0.20 


0.22 


0.00 


-0.22 


0.91 


0.90 


B^a K*^K*° 


0.19 


0.20 


0.00 


-0.22 


0.91 


0.90 



From Table I, we can see that the branching ratios for them are expected to be quite small, of 
order 10~ , since are the pure b ^ d penguin dominated decays. The subleading 

diagrams may lead to the significant CP violations in the most B decays. As 

B^ — * K^K^ decays involved only non-factorizable annihilation contributions, their branching 
ratios are much smaller than those of i? ^ K^K^, K^K^ decays, we would not study the 
B^ K^K^ modes in this paper. It should be noted that the amplitude for B^ K^K*^ is 
not simply related to that for B^ K^K*^ since the spectator quark is part of the in the 
latter decay, while in the former in the K*^. 
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Although recent experimental results in B ^ seem to be roughly consistent with 

the SM predictions, there are still windows for NP in these processes. We now turn to the RPV 
effects in B ^ }{(*)}{(*) decays. There are five RPV coupling constants contributing to the 
eight B Ki*)K^*^ decay modes. We use B, A^p and the experimental constraints shown in 
Eq. ()29|) to constrain the relevant RPV parameters. As known, data on low energy processes 
can be used to impose rather strictly constraints on many of these couplings. In FigEl we 
present the bounds on the RPV couplings. The random variation of the parameters subjecting 
to the constraints as discussed above leads to the scatter plots displayed in Fig0 




Figure 5: The allowed parameter spaces for the relevant RPV couplings constrained by -B — > 
K^*^K^*\ and (pRpv denotes the RPV weak phase. 

From Figl^l we find that every RPV weak phase has two possible bands, one band is for 
positive value of RPV weak phase, and another for negative one. We also find the magnitudes 
of the relevant RPV couplings have been up limited. The upper limits are summarized in Table 
II. For comparison, the existing bounds on these quadric coupling products IHH] are also listed. 
Our bounds on |A'j]^3A^22l5 |A'ii2Ai32l and |A^23Ai2il stronger than the existing ones. 

Using the constrained parameter spaces shown in FigEl one can predict the RPV effects 
on the other quantities which have not been measured yet in i? ^ /^(*)^(*) decays. With the 
expressions for B, A'^p and J'l at hand, we perform a scan on the input parameters and the new 
constrained RPV coupling spaces. Then the allowed ranges for B, Acp and are obtained 
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Table II: Bounds for the relevant RPV couplings hy B ^ /^(*)/^(*) decays for 
100 GeV sfermions and previous bounds are listed for comparison. 



Couplings 


Bounds [Process] 


Previous bounds [Process] 


1 \" \"* 1 


< 2.9 X 10-3 [B - 


^ _/^(*)_^(*)] 


<5.xl0-3 [B^KK] r^t 
<6.xl0-5 [B0^4>-7TO,4,<j)] |a| 


'^il3'^i22 


< 2.2 X 10-3 [B - 


^ _/^(*)_^(*)] 


< 2.9 X 10-^[B KK] m 


'^i22-^'i31 


< 1.7 X 10-3 [B - 




< 1. X 10-^ [KK]^ 


'^'il2-^i32 


< 3.0 X 10"^ [B - 




< 4. X 10"^ [5° ^ 07r°] m 


■^'123-^^1 


< 3.0 X 10-^ [B - 




< 4. X 10"^ [5° ^ 07r°] m 



Table III: The theoretical predictions for B (in unit of 10 A'^p and /l base on the RPV 
SUSY model, which are obtained by the allowed regions of the different RPV couplings. 





\ // \ //* 

'^j23'^il2 


'^il3'^i22 


^^422^^31 


'^il2'^i32 


'^i23'^i21 


BiB+ K*+K°) 


[0.0052, 7.8] 


[0.013,5.5] 


[0.0059,6.4] 


[0.056, 1.4] 


[0.064, 1.3] 


B{B+ K+K*^) 


[0.071,5.3] 


[0.056,5.3] 


[0.0096, 5.3] 








[0.0060, 7.5] 


[0.011,5.1] 


[0.0053,6.1] 


[0.049, 1.5] 


[0.054, 1.2] 


B{B^^ K^K*^) 


[0.069,5.0] 


[0.050,5.1] 


[0.0093,5.0] 






B{B+ K*+K*^) 


[0.087, 19] 


[0.041,23] 


[0.029, 16] 






B{B^ K*^K*^) 


[0.080, 17] 


[0.039, 22] 


[0.027, 15] 






A'^^piB^, - K^K^) 


[-0.75,0.57] 


[-0.19,0.44] 


[-0.18,0.47] 


[-0.18,0.47] 


[-0.18,0.50] 


■^CpiB+ K*+K°) 


[-0.19,0.17] 


[-0.32,0.17] 


[-0.42,0.47] 


[-0.99,0.99] 


[-0.98,0.76] 


A'^UK K+K*^) 


[-0.63,0.63] 


[-0.38,0.47] 


[-0.65,0.38] 








[-0.28,0.19] 


[-0.33,0.17] 


[-0.28,0.80] 


[-0.99,0.99] 


[-0.99,0.73] 




[-0.76,0.62] 


[-0.38,0.48] 


[-0.39,0.40] 






Af^'piB^ K*+K*^) 


[-0.63,0.30] 


[-0.26,0.25] 


[-0.77,0.32] 








[-0.46,0.38] 


[-0.26,0.25] 


[-0.77,0.32] 






fL{B+ ^ K*+K*'') 


[0.72,0.97] 


[0.59,0.95] 


[0.74,0.93] 






fUBj ^ K*0K*^) 


[0.72,0.97] 


[0.59,0.95] 


[0.74,0.93] 
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with five different RPV couplings, which satisfy all present experimental constraints shown in 

We obtain that the RPV effects could alter the predicted B and Aq<p significantly from their 
SM values. For decay modes, which have not been measured yet, their branching ratios can be 
changed one or two order(s) of magnitude comparing with the SM expectations, 

9. X 10"^ < B{B K+K*\ K^K*^) < 5. x 10~^ 
5. X 10"^ < B{B K*+K^, K*^K°) < 8. x 10-^ 
3. X 10"^ < B{B K*+K*^, x*°x*o-^ < 2. x 10-^ (30) 

especially, the upper limit of B{B — K*^K*^) < 2. x 10^^ which we have obtained is smaller 
than the experimental upper limit < 7. x 10^^. For Aqp, the RPV predictions on two decays 
B K*+K*^, g^rg 

A'^SUB IC+K*^) < 0.32, ^cp(^ ^ < 0.38, (31) 

and there are quite loose constraints on the direct CP asymmetries of the other five decays B — >• 
K^K^, K*+K°, K+K*^, K*°K^, K^K*^. But the RPV effects on the hiB K*+K*^, x*^x*o^ 
are found to be very small, fiiB j{*o ^ j^*o j^*o^ found to lie between 0.7 and 1, 

and these intervals are mainly due to the parameter uncertainties not the RPV effects. So 
we might come to the conclusion, the RPV SUSY predictions show that the decays B — >■ 
K*'^K*^, K*^K*^ are dominated by the longitudinal polarization, and there are not abnormal 
large transverse polarizations in Bu^ K*K* decays. The detailed numerical ranges which 
obtained by different RPV couplings are summarized in Table III. 

In FigslUHTUl we present correlations between the physical observable B, A^p, Jl and the 
parameter spaces of different RPV couplings by the three-dimensional scatter plots. The more 
information are displayed in FigslHHTIH we can see the change trends of the physical observable 
quantities with the modulus and weak phase of RPV couplings. We take the first plot in 
Figini as an example, this plot shows that B{B K*^K^) change trend with RPV coupling 
A^23'^ii2- We also give projections on three vertical planes, the \)^'l2'iKi2\'4'ijip plane display the 
allowed regions of y^i2z^ii2 which satisfy experimental data in Eq. ()29|) (the same as the first plot 
in FigEI). It's shown that B{B K*+K°) is increasing with |A.23A.;*2| on the B{B K*+K^) 
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-\K2zKi2\ plane. From the B{B K*+K^)-(j)i^ plane, we can see that B{B K*+K^) is 
increasing with 10^^!. Further refined measurements of B{B K*'^K'^) can further restrict 
the constrained space of A^23Afi*2! whereas with more narrow space of A"23A"*2 more accurate 
B{B K*+K°) can be predicted. 

The following salient features in Figs IMTOl are summarized as following. 

• Figl displays the effects of RPV coupling A^'23A-i2 on B, A'^i'p and in 5 K'^*'>K^*\ 
The constrained \^i23Ki2\~^^p plane shows the allowed range of \i23^ii2 as in the first 
plot of FigEl The six B{B K*+ K'^ , K+ K*^ , K*'^ K'^ , K*'^ , K*+ , K*'^ K*'^) have 
the similar change trends with |A'/23Afj*2l and \(f)^^ |, and they are increasing with |A'/23Afj*2l 
and |. IAqp^B K'^K^)\ are increasing with 10^^ |, but |A'/23A'/j*2| has small effect on 
A'^'piB K^K^). The two \A'}i'p{B K+ K*^ , K*^)\ tend to zero with increasing 
l'^'/23^'/i2l and |0^J. The other four \A'^SpiB K*+ , K*^ K° , K*+ K*^ , K*^ K*^)\ tend 
to zero with increasing \(f)^^ |, and they could have smaller ranges with |A'/23Afi2l- The RPV 
effects on the ^ k*+ K*'^ , K*^ K*'^) are very small, and ^ k*+K*^, K*^K*'^) 
are found to lie between 0.72 and 0.97. 

• The effects of Xii3X'i*22 on B, Acp and fi are exhibited in Fig0 The constrained |A^]^3A'j22|- 
(pj^p plane is the same as the second plot in Fig El The effects of A'ji3A'j22 on B, A^p and 
fi are similar to Af23Afi2 shown in FigEl 

• In FiglHl we plot B, A^p and fi as functions of A^22Ai3i- The constrained |A^22A'i3i|- 
(f)^p plane is the same as the third plot of FigEl The six branching ratios are in- 
creasing with |A^22'^'t3il and decreasing with \Aqp{B — > K^K'^)\ is unaffected 
by |A'j22A^3i|, but the other six direct CP asymmetries could have smaller ranges with 
|Ai22Ai3il- lAcpiK*"^ K° , K*'^ K'^)\ tends to zero with decreasing l^^^l, however, has 
small effect on Af^^B ^ K°K°, K+K*^, K'^K*^, K*+K*^, K*°K*°). The A^22A-3i effects 
on the hiB k*+ K*'^ , K*'^ K*°) are small. 

• RPV coupling A'^i2A-*32 contributes to the decays B K+K°, K^K^, K*+K^, K*^K^, and 
the effects are shown in FiglHl The constrained |A^]^2Ai32|-09p plane is the same as the 
fourth plot in Fig0 We can see that B{B K*+K^,K*^K^) are rising with |A'ii2A-32|, 
and unaffected by . Acp{B K^K^) is steady against |A-i2A-32|, and \Acp{B — > 



K*+K^,K*'^K°)\ could have smaller ranges with {X'a^^^^l- ^cp(^ ^ K^K°,K*+K^, 
K*^K^) are becoming large with increasing of 10^^ |. 

• ^i23^i2i also only contributes to the decays B K+ K'^ , K° , K*+ , K*^ , and its 
effects are shown in FigEl The constrained |A^23-^'i2il"09p plane is the same as the last 
plot in Figini B{B K*~^ , K*^ K^) are increasing with |A-23A-2i|, and unaffected by 
-^cpiB K^K^) is steady against lA'^ssA^il. and \Af.'p{B K*+ , K*^ K^)\ 
could be varied in small ranges with |A^23A-2il- ^cp(^ K°K°) is decreasing with 
109^ I, but Afi'piB K*+K^, K*^K^) are increasing with 10^^ |. 

The predictions of B and are quite uncertain in the RPV SUSY, since we just have few 
experimental measurements and many theoretical uncertainties. One must wait for the error 
bars to come down and more channels measured. With the operation of B factory experiments, 
large amounts of experimental data on hadronic B meson decays are being collected, and 
measurements of previously known observable will become more precise. From the comparison 
of our predictions in Figs IMTOl with the near future experiments, one will obtain more stringent 
bounds on the product combinations of RPV couplings. On the other hand, the RPV SUSY 
predictions of other decays will become more precise by the more stringent bounds on the RPV 
couplings. 

5 Conclusions 

In conclusions, the pure penguin B decays are very important for understanding 

the dynamics of nonleptonic two-body B decays and testing the SM. We have studied the 
B K^*)k'^*) decays with the QCDF approach in the RPV SUSY model. We have obtained 
fairly constrained parameter spaces of the RPV couplings from the present experimental data 
of B ^ I&)K'^*) decays, and some of these constraints are stronger than the existing ones. 
Furthermore, using the constrained parameter spaces, we have shown the RPV SUSY expec- 
tations for the other quantities in i? ^ /('(*)^(*) decays which have not been measured yet. 
We have found that the RPV effects could significantly alter B and from their SM values, 
but fL{B K*^K*^, ]^*0]^*o^ g^j,g T^Q^ significantly affected by the RPV effects and the decays 
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B — > K*~^K*^, ^*o^*o g^jjj dominated by the longitudinal polarization. We also have pre- 
sented correlations between the physical observable B, Af}p, /l and the constrained parameter 
spaces of RPV couplings in FigslHlfTIH which could be tested in the near future. 
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Appendix 



A. Correction functions for B M1M2 decay at ag order 

In this appendix, we present the explicit form for the correction functions appearing in the 
parameters and It's noted that in B ^ PV decays, —* ^]^{u) if M a vector 

meson. 



A.l The correction functions in B ^ PP,PV decays 

• One-loop vertex correction function is 

= 121n— - 18 + 3 f\u(^— ^ In u- in] (^M^iu). (32) 
fj, Jo \ 1 — u J 

• The hard spectator interactions are given by 

=iVc m|Fo^-^^-K,j io T^^^^Vo T^-^^^^io m7*^ H'^''^ 

if M2 is a pseudoscalar meson, and 



M Jo V 
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if M2 is a vector meson. 

Considering the off-shellness of the gluon in hard scattering kernel, it is natural to associate 



a scale Hh ~ \J AqcD^b , rather than /i ~ nih. For the logarithmically divergent integral, we will 
parameterize it as in [T7]: Xjj = Qdu/u = —ln{AQCD/T^b) + Qh^^'^" ^b/^QCD with {qh^'Ph) 
related to the contributions from hard spectator scattering. In the numerical analysis, we take 
Aqcd = 0.5GeV, {gt, (pu) = (0, 0) as our default values. The same as in i? ^ VV decay. 
• The penguin contributions at the twist-2 are described by the functions 



pp 



CiGm,{sp) + Cs GmM + + (C4 + Ce) 



[nf-2)GMM 



2n/ 



U 



PIhJ = {G, + NcG2)Gm,{s,)- G!;^^^ J\u^ 



u 



(35) 



where Uf = 5 is the number of quark flavors, and s„ = 0, Sc = (mc/mfe)^ are mass ratios involved 
in the evaluation of the penguin diagrams. The function Gm2{s) is defined as 



nib 



Gm2{s] 

• The twist-3 terms from the penguin diagrams are given by 



— I — In 1-4 du dx xxln{s — xxu — ie)^M2\'^)■ 
3 3 ft Jo Jo 



(36) 



Pk,3 = C,Gm2{sp) + Gs Gm2{0) + Gm2{1) + (C4 + C6) (n^-2)GM,(0) 



GM,is,)+G 



M2 



2n, 



Pf}'^^ — {Gi + NcG2)Gm2{sp) — "iGjl^^, 



M2,3 



(37) 



with 



— \ — In h4 du dx xxln{s — xxu ~ ie)^^^'^{u) 

3 3 " 



(38) 



jj. Jo Jo 

if M2 is a pseudoscalar meson, and we omit the twist-3 terms from the penguin diagrams when 
M2 is a vector meson. 

• The weak annihilation contributions are given by 



Al(Mi,M2) ^ ^^2(^1,^2) ~ vra. 



18 - 4 



TT 



A^(Mi, M2) ^ QTias{r'^' - <^^) - 2Xa + yj , 
A{(Mi, M2) ^ 67ra,(rf + ) (2X1 - 
A((Mi,M2) = 0, 4(Mi,M2) = 0. 



(39) 
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when both final state mesons are pseudoscalar, whereas 

Ai(Mi, M2) ^ -v4^(Mi, M2) ^ ISvra, f - 4 + 



TT 



2' 



A'3(Mi, M2) ^ 67ra,r^^^^ (^X^ - 2Xa + yj , 

4(Mi, M2) ^ -evra.rf (2X1 " Xa) , 

M2) = 0, M2) = 0. (40) 

when Ml is a vector meson and M2 is a pseudoscalar. For the opposite case of a pseudoscalar 
Ml and a vector M2, one exchanges r^^ ^ r^^ in the previous equations and changes the sign 
of 4. 

Here the superscripts i and / refer to gluon emission from the initial and final state quarks, 
respectively. The subscript k of A^jl^ refers to one of the three possible Dirac structures Fi ® F2, 
namely k = I loi {V - A)®{V - A), k = 2 for {V -A)®{V + A), and A; = 3 for {-2){S-P)® 
(S + P). Xa = Jq dn/n is a logarithmically divergent integral, and will be phenomenologically 
parameterized in the calculation as Xh- As for the hard spectator terms, we will evaluate the 



various quantities in Eqs. (jHIIjl and at the scale fih = y-^QCDfnt,. 
A. 2 B ^VV decays 

In the rest frame of B system, since the B meson has spin zero, two vectors have the same 
helicity therefore three polarization states are possible, one longitudinal (L) and two transverse, 
corresponding to helicities A = and A = ± ( here Ai = A2 = A). We assume the Mi(M2) meson 
fiying in the minus(plus) z-direction carrying the momentum pi(p2), Using the sign convention 

g0123 ^ Yiaye 



A 



M1M2 



ho, 



(41) 



where for A = and h± for A = ± 

nih , „ , \ 1 ~ 2m 



• V^A (±1) contain the contributions from the vertex corrections, and given by 



1/° (a) = 121n— - 18+ / du^fHu) [3 Inw-SiTr , (42) 

H JO " V 1 — M 



l/4(a) = 121n^-18+ (" du U^'Hu) fsl^ In n - 3zvr^ 

jji Jo \ A j \ 1 — u J 
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• The hard spectator scattering contributions, exphcit calculations for H^j^j^j^{a) yield 



Nc ho Jo ^ JO V Jo u 



47r^ 2ifBfMjM2mM, 
Nc mBh± 
1 



Jo t Jo 



dv- 



y2 



l(a)M2 I \ 

X / du\^^''^u)-'^^^ 







X / dt^dn 
Jo 



(v)M: 



4 

/(a)Mi 



e Jo 

4:71^ 2ifB /m, /m2 rriM, ^^^liO 
Nc m%h± Jo ^ 



V) \ I (v)M2/ \ , ^' 



l{a)M2 



{u) \ U + V 



(43) 



with V = 1 — V, when the asymptotical form for the vector meson LCDAs adopted, there will 
be infrared divergences in i?MiA4'2- OEHl, we introduce a cutoff of order kqcolmh and 

take IS.QCD = 0.5 GeV as our default value. 

• The contributions of the QCD penguin-type diagrams can be described by the functions 



^M2,2 



pX,p,EW 
^M2,2 



q—u 



+\Co[e,Glj^{s,) + e^G\,^{s,)\ + ^(Cg + G^o) E e<?' 



2 '-'77 '-^g' 



(44) 



G^2(«) = ^ + ^ln^ + 4 /'dn^f (u)g(u,.). 



2 2, mfc /-i 
- + -ln^ + 2 / dM 

O O /i JO 



l{a)M2 



U 



with the function g(u, s) defined as 

g(u, s) 



dx xa; In is — xxu — ie) 



{u,s), 



(45) 



(46) 



We omit the twist-3 terms from the penguin diagrams for B VV decays. 

• We have also taken into account the contributions of the dipole operator Og^, which are 
described by the functions 

1 2^1'' (u) 







- du 




Jo 




/■I du 




Jo u 



1-u 



-ugy [u) T 



UK ' Uu) 



+ dv[^f'{v) 



{y)M2 



(a)M2 



U 



(47) 
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here we consider the higher-twist effects /c'* = uEnt + kj_ + -^n'^ in the projector of the vector 
meson. The = in Eq. (fTTjl j2Sl IHO] if considering the Wandzura- Wilczek-type relations 

I2n!. 

We have not onsidered the annihilation contributions in B VV decays. 
A. 3 The contributions of new operators in RPV SUSY 



Compared with the operators in the 'H^^, there are new operators {q2(l3)v±A{bqi)v+A in the 
• For B PP, PV decays, since 



{P\ gi7/.(l - 75)^2! 0) = -{P\ 517^(1 + 75)92! 0) = -{P\ gi7/.75g2| 0), 
(P| g7^(l - 75)6 I B) = {P\ g7^(l + 75)6 | B) = {P\ q^,b \ B), 
{V\ qil,{l-lM 0) = {V\ qil,{l+lM 0) = {V\ gi7^g2| 0), 
{V\ q^,{l - 75)6 I B) = -{V\ 57^(1 + 75)& I B) = -{V\ q^,^,b \ B), (48) 

the RPV contribution to the decay amplitude will modify the SM amplitude by an overall 
relation. 

• For B — > VV, we will use the prime on the quantities stands for the {q2q3)v±A(bqi)v+A 
current contribution. In the NF approach, the factorizable amplitude can be expressed as 

^M.M, = (M2|fe7M(l-«75)g3)|0)(Mi|(67^(l + 75)gi)|5). (49) 
Taking the Mi(M2) meson flying in the minus(plus) z-direction and using the sign convention 

g0123 ^ _1 \]_Q^^Q 



A' 



(ml - ml,^ - minims + m^JAf^^^^ {ml, J - ^^^^A^^'^^ {ml, J 



(50) 



The vertex corrections V^}^{a) and the hard spectator scattering corrections H'^^j^,j^{a) as 
follows: 

v't^iAa) = -121n— + 18-6(l + a)- /'du$f^(u) fs^^^lnu-SiTrV 

Vl±(a) = -12 In ^ + 18 - 6(1 + a) - du Ul>''\u) ± "g^"^^""^ (") (^3!^ \nu- 3..) , 
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Nc K Jo i Jo V Jq 

""-^-^^"^ = ^'^'Vo '^^io '"^^ 



A ( (v)M2( \ ag[^ {u)\ Att^ 2tfBfMifM2mMimM2 f ^c^iiO 
dulgy {u) + +— / — ^ 

X / d.d. g^l^^'^iv) t'-^ ^ ) g'l'''^{u) ±^ ^ 1 



B . The amplitudes in the SM 



K°K*° 



04 - -aio + r 
2 ^ 



(ae - -as) 



V2 



tb^td 



1 K°, 1 N 

a4--a.io + r^ [ag - -as) 



} 



-Vtiytd 



V2\ 



tb^td 



b3{K", K") + b4{K", K") + b4(K", K") 



-^bf"iK'^,K°) - ^bf^(K°,K°) - ^bl"' {K° , K'^) 



a4 - -aio - r. 



(ae - -as) 



^fBfKfK' {v:,V^db2iK*+,K0)-V,lVta [b3{K*+ , K°) + b^-^ {K*+ , RO)]} 



tb^td 



14, - -aio 



} 



V2 

^fsfKfK' {v:,V„ab2{K+,K''') - v:,Vu [63(i^+,i?*°) + f)r(/f+, A-*")] } : 
V2 



{-Vtl 



tb^td 



a4 - -aio - ("6 - -as) 



^fBfKfK' {- 



bi{K'°, K°) + b4{K*°,K°) + b4{K°, K'°) 
1 . 



.0 ^?0^ 



2 

G F 

7f 



} 



^tb^td 



04 - -aio 



b-i{K°,K*'^) + b4{K°,K*°) + b4{K*°, K°) 



■hl'^{K°,K*°) - if)|"'(A'°,_ft:*0) - ^bl'^ {K*° , K°) I 



tb^td 



a4 - -aio 



a4 — 2'^10 



Here we have not considered the annihilation contributions in i? — VV decays. 



C . The ampHtudes for RPV 

-■■ ( Ml Ml* 

A'^ p {B+ ^ K+ K°) = 
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8m? 
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-4/00 ^ 



'^il2^»32 
8m? 



^i23^i21 
8m? 

•^Li 



^i23^il2 „-4/0o P - _i_ 



16m? 



_i23_il2_ 
16m?. 



-4/00 p 



8m? 

\ / \/* 

^il3 i22 

8m? 

I^Li 

-8/00^^0 
X 



8m? 

\/ \ /* 
^113 t22 

8m? 



^i22^i31 
8m? 

I^Li 



Sm'i 

•^Li 

\f \ 

i22^i31 

8m? 



-8/00 2, 



K+K'O 



^-8/0O(_i^.„^„) 



-8/00 r' /l' 
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In the ^^J', F^j^j^i^ and Lm^^.j^ are defined as 
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LmiM2 

PP, PV decays, and 



^ ]_ as Cf 

Nc 47r A^c 

A^c I 47r Nc 



VkPi + HmiM2 

12 + Va/2 + HmiM2 



F' 



Tl 
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1 - 
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AtiNc 
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Us Cf 




{- 


AtiNc 


1 




as Cf 




{- 


A-K Nc 



ylX 



i2 + v^Xi^) + H''m 



i2 + v^4( 



-l)+^^M.A/2(-l) 



B — > VV decays. 
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Figure 6: The effects of RPV coupling AfsgAffg in B ^ k^*)k^*) decays. 
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Figure 7: The effects of RPV coupling A'^^gAfsa in 5 ^ k'^*)k'^*) decays. 
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Figure 8: The effects of RPV coupling ^22^^,1 in -B ^ k'^*)k'^*) decays. 
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Figure 9: The effects of RPV coupling A'^igAfga mB ^ k^*)k'^*) decays. 




31 



